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Since its inception, the quark model has had great suc- 
cess in describing the spectroscopy of hadrons. The quark 
model has been successful for the B mesons, where all 
of the ground states have been observed The spec- 
troscopy of c baryons also agrees well with the quark 
model, and a rich spectrum of baryons containing b 
quarks is predicted [2j. However, direct observation of b 
baryons has been limited to a single state, the A& (quark 
content udb) until recently. Evidence for b baryons 
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that also contain a strange quark was shown from LEP [3] 
through partial reconstruction of decays containing elec- 
trons and muons. Recent results from the Tevatron on 
the Tit, states (quark content uub, ddb) 0] and Hr (quark 
content dsb) [5J are beginning to subject the b baryons 
to closer examination. 

In this Letter, we report the observation of a heavy 
baryon and measurement of its mass. The decay proper- 
ties of this state are consistent with the weak decay of a b 
baryon, and we interpret our result as the observation of 
the S7~ baryon. This observation is made in pp collisions 
at a center of mass energy of 1.96 TeV using the Collider 
Detector at Fermilab (CDF II), through the decay chain 
— » J/ipa~, where J/tp — » S~ — > Air~, and 

A — » pir~ . Charge conjugate modes are included implic- 
itly. This measurement is based on a data sample with 
an integrated luminosity of 1.9 fb _1 . 

The CDF II detector has been described in detail 
elsewhere This analysis primarily relies upon the 
tracking and muon identification systems. The track- 
ing system consists of a silicon microstrip detector and 
an open-cell drift chamber (COT) that operate inside 
a solenoid with a 1.4 T magnetic field. Muon candi- 
dates from the decay J/tp — > fi + /jT arc identified by two 
sets of drift chambers located radially outside the electro- 
magnetic and hadronic calorimeters. The central muon 
chambers cover the pseudorapidity region \t)\ < 0.6 Q, 
and are sensitive to muons with transverse momentum 
Pt > 1.4 GeV/c. A second muon system covers the re- 
gion 0.6 < \r)\ < 1.0 and detects muons having px > 2.0 
GeV/c. Muon triggering and identification are based on 
matching tracks measured in the muon system to COT 
tracks. 

The analysis of the data begins with a selection of well- 
measured J/tp — > candidates. The trigger require- 
ments are confirmed by selecting events that contain two 
oppositely charged muon candidates, each with matching 
COT and muon chamber tracks. We also require that 
both muon tracks have associated measurements in at 
least three layers of the silicon detector and a two-track 
invariant mass within 80 MeV/c 2 of the world-average 
J /if) mass Tj. This data sample provides approximately 
15 million events containing J/tp candidates, measured 
with an average mass resolution of 20 MeV/c 2 . 

The reconstruction of S~ candidates uses all additional 
tracks found in each selected J/tp event. Pairs of oppo- 
sitely charged tracks are used to identify A decay candi- 
dates. The proton (pion) mass is assigned to the track 
with the higher (lower) momentum. This mass assign- 
ment is always correct for A — > p tx~ candidates used in 
this analysis because of the kinematics of A decay and the 
lower limit of «200 MeV/c in the transverse momentum 
acceptance of the tracking system. The A mass is mea- 
sured with a resolution of 2.5 MeV/c 2 . All intersecting 
pairs of tracks with an invariant mass within 10 MeV/c 2 
of the world average A mass [l| have their track parame- 
ters recalculated according to a fit where the momenta of 
the two tracks are constrained to the A mass. The decay 



vertex is used to calculate the A displacement from the 
beam-line in the direction of the track pair's transverse 
momentum. The background due to tracks originating 
from the primary vertex is reduced by requiring this dis- 
placement to exceed 1.0 cm. For candidates that satisfy 
these requirements, the remaining tracks are assigned the 
pion mass, and A 7r~ combinations are identified that are 
consistent with the decay process S~ — > A-7r~. In order 
to obtain the best possible A ir~ mass resolution, the re- 
construction uses a fit on the three tracks that simultane- 
ously constrains the A decay products to the A mass, and 
the A trajectory to intersect with the helix of the S~ de- 
cay pion. For all S~ candidates, the reconstructed decay 
position of the A candidate is required to be radially dis- 
placed at least 1.0 cm with respect to the reconstructed 
decay vertex of the S~ candidate. 

The majority of S~ candidates have pr > 1.5 GeV/c. 
This, along with the long lifetime of the S~ (cr = 4.9 
cm) [ljj, results in a significant fraction of the S~ can- 
didates having decay vertices located several centimeters 
radially outward from the beam-line. Therefore, we are 
able to refine the S - reconstruction by making use of the 
improved determination of the trajectory that can be ob- 
tained by tracking the S~ in the silicon detector. The S~ 
candidates have an additional fit performed on the three 
tracks that simultaneously constrains both the A and S~ 
masses of the appropriate track combinations, and pro- 
vides the best possible estimate of the S~ momentum 
and decay position. The result of this fit is used to de- 
fine a helix that serves as the seed for an algorithm that 
searches for silicon detector hits associated with the 5~ 
track. We retain for further analysis all S~ candidates 
with measurements in at least two layers of the silicon 
detector. This technique provides excellent impact pa- 
rameter resolution for the S - track (average of 60 /an), 
and has been used previously Q . The A ir~ invariant 
mass spectrum of all combinations that satisfy these re- 
quirements is shown in Fig. [T] Approximately 23,500 
S~ candidates above the combinatorial background are 
identified. The search includes the subset of these 
combinations with an invariant mass within 10 MeV/c 2 
of the world average S~ mass [l|. 

The mass resolution for the J/tp S~ final state is stud- 
ied with a Monte Carlo simulation that generates b 
quarks according to a next-to-leading-order calculation 
and produces events by simulating b quark frag- 
mentation [nj. The decay -> J/tpE~ is simulated 
with EvtGen [ll[. The generated events are used as 
input to the detector and trigger simulations based on 
a GEANT3 description [12j and processed through the 
same reconstruction and analysis algorithms used for the 
data. Analysis of the simulated events shows that 
a 10% improvement in mass resolution can be obtained 
if the momenta of the S~ decay products are allowed to 
vary in the fit of the candidate, rather than simply us- 
ing the S - track. Consequently, a procedure that simul- 
taneously fits the five tracks of the final state, constrains 
the three vertices of the decay chain to the appropriate 
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topology, and constrains the masses of the J/ip, E - , and 
A to their world average masses, is used to provide the 
best estimate of the J/ip S~ mass. The average S,J~ mass 
resolution obtained from simulated events is found to be 
approximately 15 MeV/c 2 . In particular, we note that 
the J/ip S~ invariant mass resolution is comparable to 
the mass resolution obtained with the CDF II detector 
for other B hadrons with a J/ip in the final state |l3|. 

The selection used to isolate the — > J/ip S~ decay 
process is guided by the properties of other B hadrons 
that include a J/ip in the final state. The important 
ones include the lifetime of the ground-state B hadrons 
and the energy available in the decay. The B ± ,B°, and 
B a mesons and Af, baryon, all appear to have lifetimes 
dominated by the weak decay of the b quark. We expect 
the same to hold true for the and for its lifetime to 
be comparable to these states. In addition, these parti- 
cles all decay to final states J/ipX, where A" is a single 
hadron. In these decays the momentum carried by the 
decay products in the rest frame of the B hadron falls in 
the fairly narrow range of 1570— 1744 MeV/c 2 for decays 
where only the lightest decay products are considered. 
We expect the energy released in the decay of the to 
be comparable. This expectation is also consistent with 
the range of theoretical predictions (5788-5812 MeV/c 2 ) 
for the mass [2] . Consequently, we have chosen a SjJ~ 
search strategy that will provide an optimal sensitivity 
for a final state that shares similar properties with these 
well-established B hadrons. 

We have developed the event selection by studying 
B ± — » J/ipK ± decays. This final state is identified by 
assigning the AT ± mass to all tracks not used in the J/tp 
reconstruction. Each three-track combination must sat- 
isfy a fit where the tracks are required to originate from a 
common vertex and the invariant mass of the muon pair 
is constrained to the world average J/ip mass Ap- 
proximately 30,000 S ± candidates are identified in this 
sample. Several characteristics of the final state are used 
as selection requirements to obtain a B signal with very 
little background. Minimum transverse momentum re- 
quirements on the and B candidates are used to 
suppress backgrounds from the event that are not re- 
lated to the B^ decay. The trajectory of the AT ± is re- 
quired to originate from the B decay vertex by placing 
a requirement on its impact parameter dsv(K) and as- 
sociated uncertainty <Jd sv {K) with respect to the vertex 
found in the J/ip fit. Similar impact parameter quantities 
dpv (K) and <Jd PV (K) measured with respect to the pri- 
mary vertex are used to remove tracks that originate from 
the prompt background. Reasonable vertex quality is as- 
sured by placing a minimum value on the accepted prob- 
ability P(x 2 ) of the mass- and vertex-constrained fit used 
to obtain the B ± candidate. We suppress the promptly- 
produced combinatorial background by rejecting candi- 
dates with low proper decay time, t = Lt -pr(B) r^ffgjp , 

where M(B) is the mass of the B ± candidate, pr{B) 
is the transverse momentum of the B^ candidate, and 
Lt is the transverse displacement of the B decay ver- 



tex from the beam-line. A requirement on proper decay 
time uncertainty a t removes poorly-reconstructed com- 
binations. We also reject combinations that are inconsis- 
tent with having originated from the beam-line by requir- 
ing a small magnitude of the impact of the B candidate, 
dpv(B) = Lt x Pt(B)/\pt{B)\, and a small angle (3 be- 
tween Lt and pr(B). 

This analysis uses a two-step selection procedure. Fi- 
nal selection criteria are listed in TableU We first impose 
the "standard" selection requirements listed there and 
retain all J/ip K combinations that satisfy them. Any 
combination that fails only one of the "standard" selec- 
tion requirements is also allowed into the final sample 
if it satisfies both of the pt requirements of the "high- 
Pt" selection requirements and fails no more than one of 
the other requirements in this set. The combination of 
"standard" and "high-p^" requirements reduce the back- 
ground in the B^ — > J/ip sample to approximately 
400 combinations, while retaining a signal of 16,000 B ± 
candidates. 

As is done for the reconstruction of the B ± , the treat- 
ment of the J/ipS~ candidates requires a mass- and 
vertex-constrained fit on the muon candidates and the 
S~ track. The selection criteria in Table U are applied to 
the J / ip S~ sample where we simply exchange K~ for S~ 
and B ± for where appropriate. Combinations that 
satisfy these requirements form the final set of candi- 
dates. The invariant mass of each candidate is obtained 
with the full five-track fit, and the resulting J/ip S _ mass 
distribution is shown in Fig. O 



TABLE I: Selection variables and requirements for the "stan- 
dard" selection and "high-pr" selection as described in the 
text. Here l 'K" refers to the third track combined with the 
J/ip and is a K or E~ candidate for the B or Hr candi- 
dates, respectively. Similarly, "B" refers to either B or , 
as is appropriate. 



Selection variable 


Standard 


High-pr 


p T (K) 


> 1.7 GeV/c 


> 2.5 GeV/c 


Pt(B) 


> 5 GeV/c 


> 6 GeV/c 



\dsv(K)\ 
\dpv(K)\/vd PV (K) 

P(x 2 ) 

\d PV {B)\ 
ct 
cat 





< 100 /an 
> 2.5 

> 0.1% 

< 75 nm. 
> 80 fim 

< 30 /zm 

< 0.4 radians 



> 3 
> 1% 

< 60 /im 
> 100 ^m 

< 25 fxm 

< 0.3 radians 



The mass resolution estimate for the 5 fe implies that 
more than 95% of a signal will occupy an invariant 
mass bin with a width of 75 MeV/c 2 . The data shown in 
Fig. O contain 18 candidates in the 75 MeV/c 2 range of 
5750 - 5825 MeV/c 2 . We model the combinatorial back- 
ground by considering candidates in the range 5700—6500 
MeV/c 2 ; the data yield 23 candidates in this range. No 
events contribute multiple candidates. The upper limit 
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of this range is chosen arbitrarily, and has no impact on 
the result. The lower limit is chosen to avoid partially 
reconstructed H^~'° — > Jj^'Er X decays, where X rep- 
resents additional undetected particles. We assume that 
the mass distribution of the combinatorial background is 
uniform and that the occupancy due to background com- 
binations in any particular 75 MeV/c 2 mass bin within 
the 800 MeV/c 2 search range can be described by a bi- 
nomial distribution, with a single event probability given 
by the ratio of the two mass ranges used. The prob- 
ability that the number of candidates observed in the 
5750 — 5825 MeV/c 2 mass range is due to a background 
fluctuation is estimated as the binomial probability of 18 
or more events from a sample of 23 total occurences and 
a single event probability of 75/800. This probability is 
6.6 x 10~ 15 , equivalent to a 7.7c variation from a Gaus- 
sian distribution. Consequently, we interpret the data 
distribution shown in Fig. [2] to be the observation of 
a resonance, with a width consistent with the detector 
resolution. Comparable distributions of J/ip A tt + and 
J/^A7r~ where the Att~ does not form a 5~, yield no 
significant enhancement at any mass within the range of 
this analysis. 

The masses and their uncertainties obtained from the 
five-track final state fit method are used in an unbinned 
likelihood fit to measure the mass. The negative log- 
likelihood function that is minimized has the form 

N 

£ = In [/GK, m , s m a?) + (1 - f)C] , (1) 

i=l 

where mj is the mass obtained for a single candidate, cr™ 
is the uncertainty on that mass as estimated from the 
track parameters for the candidate, G(rrii, too, s m a™) is 
a Gaussian distribution with average too and character- 
istic width s m cr™, and C is a constant background term. 
The quantities obtained from the fitting procedure in- 
clude /, the fraction of the events found in the signal, 
too, the best average mass, and s m , a scale factor on the 
mass uncertainty which accounts for a possible shift in 
the mass uncertainty. This procedure yields a best es- 
timate for the Z~ mass of 5792.9 ± 2.5 MeV/c 2 . The 
uncertainty scale factor is determined to be 2.0 ± 0.4, 
which is consistent with the value of 1.6 ± 0.2 obtained 
with simulated candidates. The signal fraction is cal- 
culated to be 0.76 ± 0.09, giving a yield of 17.5 ± 4.3 S~ 
candidates in this data sample. The projection of the fit 
function is superimposed on Fig. [2] 

Several systematic effects have been considered for 
their impact on the mass measurement. The over- 
all momentum scale of the tracking system is established 
by calibrating with the well-measured J ftp, ip(2S), and T 



states [13]. The momentum scale uncertainty contributes 
a systematic uncertainty of ±0.4 MeV/c 2 on the !E^~ 
mass measurement. Alignment and material distribution 
of the tracking system contributes an additional ±0.6 
MeV/c 2 of uncertainty. The method requires knowledge 
of the mass of the final state S~ , a quantity that is known 
to ±0.13 MeV/c 2 [l|. The mass uncertainty of the J/tp 
is included in the momentum scale calibration, and the 
mass uncertainty of the A (±0.006 MeV/c 2 ) is negligible 
for this analysis. Finally, the largest systematic variation 
seen on the mass measurement occurs when alternative 
fitting models are used for the mass calculation. The sig- 
nal distribution has been fit with a single Gaussian where 
its parameters are allowed to vary in the fit. A double 
Gaussian is also used, where the widths are fixed to val- 
ues obtained in simulation, and only the average mass is 
allowed to vary. These different models for the probabil- 
ity distribution function of the signal create variations of 
±1.5 MeV/c 2 on the mass result. The individual system- 
atic uncertainties are combined in quadrature to obtain 
an overall systematic uncertainty of ±1.7 MeV/c 2 on the 
mass measurement. 

In conclusion, we use the CDF II detector at the Teva- 
tron to observe the in pp collisions. A signal with 
17.5 ± 4.3 candidates and a significance of 7.7er is 
seen in the decay channel — > J/ipE^. The mass of 
this baryon is measured to be 5792.9 ± 2.5 (stat.) ±1.7 
(syst.) MeV/c 2 , which is consistent with theoretical ex- 
pectations 0. The mass measurement presented is also 
consistent with the only other direct observation of this 
state Q , and represents a significant improvement in pre- 
cision. 
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FIG. 1: The invariant mass distribution of An~ combinations 
having an associated track in the silicon detector in events 
containing J/ip candidates. 
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FIG. 2: The J/xjj H~ invariant mass distribution for combina- 
tions that satisfy the selection requirements. The projection 
of the fit function is overlaid on the data. 



